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Multilayer approximation for a confined fiuid in a slit pore 
G. J. Zarragoicoechea,11'2^ A. G. Meyra,^ V. A. Kn^ 

A simple Lennard-Jones fluid conflned in a slit nanopore with hard walls is studied on 
the basis of a multilayer structured model. Each layer is homogeneous and parallel to 
the walls of the pore. The Helmholtz energy of this system is constructed following van 
der Waals-like approximations, with the advantage that the model geometry permits to 
obtain analytical expressions for the integrals involved. Being the multilayer system in 
thermodynamic equilibrium, a system of non-linear equations is obtained for the densities 
and widths of the layers. A numerical solution of the equations gives the density profile 
and the longitudinal pressures. The results are compared with Monte Carlo simulations 
and with experimental data for Nitrogen, showing very good agreement. 



I. Introduction 

The effects on phase transition of confined fluids 
in a slit-like pore have been studied by simulation 
and different theories [1-11]. In a previous work, 
we constructed a generalized van der Waals equa- 
tion for a fluid confined in a nanopore [12, 13]. 
The shift of the critical parameters was in good 
agreement with lattice model and numerical sim- 
ulation results, and the predicted critical temper- 
ature remarkably reproduced the experiment. In 
that work, we concluded that the confined van der 
Waals fluid theory seemed to work better than the 
bulk one, maybe due to the fact that the higher 
virial contributions not considered in both theo- 
ries were less important in the conflned fluid than 
in the bulk. A similar treatment was used previ- 
ously by Schoen and Diestler [14]. Following that 
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line of reasoning, here we study a simple fluid con- 
fined between two infinite parallel hard walls (slit 
pore). The walls are at a distance L apart. To 
study the conflned fiuid, we propose a multilayer 
model [15]: the fluid is distributed in n thin lay- 
ers, one beside the other. Each layer has a uniform 
density, and can be observed as a non-autonomous 
phase. A particle in a given layer interacts with 
its neighbors inside the layer, and with every par- 
ticle in the other layers. Defay and Prigogine and 
Murakami et al. have shown that, in a liquid gas 
interface, the deviation from the Gibbs' adsorption 
equation becomes practically negligible in the case 
of a two layer model [16], and that as the number 
of transition layers grows, the multilayer model be- 
comes perfectly consistent with the Gibbs' equation 
[17]. The van der Waals-like approximations made 
in developing this multilayer model theory limit its 
validity to the low density regime. 

II. Theory 

The model system consists of a fluid of Lennard- 
Jones particles confined in a slit nanopore. The 
hard walls of the pore, separated at a distance L (in 
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the X direction), have a surface area S [S — > oo). 
We divided the fluid into n layers, each layer being 
parallel to the pore walls. The layer i has Ni par- 
ticles {N = X;r=i a width L„ (L= YJt^i Lxi), 
and a volume Vi = SL^i. Then the Helmholtz en- 
ergy [18] can be written as 



A = -fcTln 



n 



(1) 



The configuration integral Zpj for a pair potential 
Vij may be approximated as 



yne--^/'=v«^y/,,dr"'(2) 



i=l 

h, = e--^/'=^-l, 

and further expanded in function of two particle 
integrals 
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The first term in Eq. ([3| stands for particles in the 
layer i. The second term comes from the interac- 
tion of one particle in layer i with one particle in 
layer j. In a compact form, and assuming that a 
layer sees three nearest neighbor layers, 



Z 



N 



(4) 
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The integrals li and 7^, for the slit pore geome- 
try and after low density approximations, can be 
analytically solved to give 

Ii = JJ /i2dridr2 w - // dridr2 
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(7) 
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In the above expressions Vy was taken to be the 
Lennard- Jones pair interaction, being e and a the 
potential parameters. The integrals li^ i+2 and 
/i. i+3 do not contain the excluded volume term be- 
cause we suppose that the layer widths are Lxi ^ ^• 
The expressions for A and B in the preceding equa- 
tions, functions of L^^^ are given in the Appendix 
A. 

The Helmholtz energy, Eq. (1) together with Eq. 
(4), has the final expression 
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The pressure tensor [12. 13] and chemical poten- 
tials are obtained from the following equations 




If the system is in mechanical and chemical equi- 
librium, the XX components of the pressure tensor 
and the chemical potentials for each layer must be 
equal. Prom these equations, giving as input the 
wall separation L and the mean density p* = pcr'^, it 
is constructed a system of (n — 1) non-linear equa- 
tions with (n — 1) unknowns (layer densities and 
widths) to be numerically solved. The low com- 
putational cost is taken for granted given that the 
code is easily written and the calculations are car- 
ried out on a Pentium 4 processor running at 2.66 
GHz. At a temperature T*=kT /e=l, we have ex- 
plored the cases with L=10a and L=15a, at dif- 
ferent mean densities. We have also compared the 
theoretical results with experimental data coming 
from studies of nitrogen adsorption in graphite slit 
pores at room temperature [19]. 

III. JMonte Ccirlo simulation 

For numerical simulations, N Lennard Jones par- 
ticles are confined between hard walls separated at 
a distance L. The unit cell is biiild up taken the 
walls to be of size Ly and in the y and z direc- 
tions respectively, directions on which the period- 
ical boundary conditions are applied. The density 
profiles and pressures were obtained taking average 
values in fluid slabs parallel to the walls. The pres- 
sure tensor was used in the simple virial form, as 
indicated in references [20, 21]. 

With T*=l, and for both slit pore widths L=10a 
and L=15(T, the size of the imit cell was set to 
Ly=Lz=30a, taking the number of particles N to 
correspond with the mean density. The range of the 
Lennard- Jones interactions was considered with a 
cutoff radius of 5a. 
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Figure 1: Density profiles for a n=9 layer model 
of a confined fluid in a slit pore (solid symbols). 
The temperature is r*=1.0 and the wall separation 
is L=10a, with mean densities p*=l/20 (circles), 
1/15 (squares), 1/10 (triangles), 1/5 (diamonds), 
and I (stars). Open symbols represent the Monte 
Carlo simulations. 




X /a 



Figure 2: zz pressure tensor component. Captions 
as in Fig. 1. 

IV. Results 

In Figs. 1 and 2, the density profiles and zz 
components of the pressure tensor are shown for 
T*=l and L=W(t. The mean densities studied are 
p*=l/20, 1/15, 1/10, 1/5, and 1/4. The agreement 
of the theoretical density profiles with the Monte 
Carlo simulations is very good. For the pressure 
there is a rather good correspondence for low den- 
sities, up to p*=l/10. For the higher densities, dif- 
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Figure 3: Density profiles for a n=13 layer niodel 
of a confined fluid in a slit pore (solid symbols). 
The temperature is T* = 1.0 and the wall separa- 
tion is L—lba^ with mean densities p*=l/10 (tri- 
angles), and 1/5 (circles). Open symbols represent 
the Monte Carlo simulations. 



ferences appear, though the tendencies are similar. 
The discrepancies come first from the low density 
approximations done to get the Helmholtz energy. 
But, while in the simulation slab particles fluctu- 
ate and at higher densities some clusterization oc- 
curs, in the theory each layer is supposed to have a 
homogeneous density which makes it hard for the 
theoretical pressures to follow those obtained by 
simulation. For the density profiles, averaging the 
number of particles in each slab evidently compen- 
sates the clusterization, and the theory gives good 
results, at least for the rather low densities studied. 
The same picture applies to the behavior of the sys- 
tem for T*=l and L=15(t, at mean densities p* — 
1/10, and 1/5, represented in Figs. 3 and 4. 

The results, as expected for hard repulsive walls, 
show a low density region next to the walls and 
an increasing density profile, with a maximum at 
the center of the slit pore. This behavior is also 
shown with density functional theory [1] and in 
other Monte Carlo simulations [2]. 

Finally, the good agreement of the theory with 
the experiment can be seen in the results shown 
in Fig. 5. In this figure, the excess number of 
molecules per unit area of pore surface F is plotted 
in function of the external pressure, at T*=3.18 
and L=4(T. These parameters approximate the ex- 
perimental values [19] r=303 K and L=1.45 nm, if 



Figure 4: zz pressure tensor component. Captions 
as in Fig. 3. 

e/k = 95.2 K and a = 3.75 A are used to character- 
ize the nitrogen. In this case, due to the size of the 
sample, n=3 layers have been used for calculation. 
F is defined as 



S 



iP* Pi) 



L 



(11) 



where Ng/p*g is the number/density of particles 
which would occupy the slit pore in the absence of 
the adsorption forces, p* and the external pressure 
are determined equating the chemical potential in- 
side the slit pore (Eq. 10 1 to the chemical poten- 



tial coming from the bulk van der Waals equation 
at the same temperature. The theoretical results 
presented here are similar to the numerical simu- 
lation results obtained by the same authors who 
have done the experiment [19]. They assume that 
the differences at higher pressures could be a con- 
sequence of the uncertainty in the determination of 
the pore geometry. 

V. Conclusions 

The application of a simple theory, with van der 
Waals-like approximations to the Helmholtz en- 
ergy, to a particular model of spatial distribution 
makes it possible to obtain analytical expressions 
for the thermodynamic quantities. The study of a 
confined fiuid in a slit pore geometry with a mul- 
tilayer approximation produces good results when 
compared with Monte Carlo simulations at low den- 
sities. The agreement with a particular experi- 
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P (Mpa) 



b = l7r;Bi = lj^^;Ai = ai + + ft + ft 



ai = -^7r;a2 = §7r;a3 



-|7r;a4 = ^tt (Al) 



A correction has been made to get good critical 
parameters for the bulk (L oo). For Argon oi = 
-5.7538 and 6=1.3538, and for nitrogen ai= -1.5955 
and 6=1.0349. 
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Figure 5: Excess number of molecules per unit area 
of pore surface F as function of the external pres- 
sure. The full line represents the experiment (digi- 
talized from Ref. [19]), and the dots are our theo- 
retical results. 
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ment on nitrogen confined in a graphite slit pore 
is remarkable, even though an excess quantity is in 
study. It may be concluded that the confinement 
reduces the importance that higher virial contribu- 
tions have on the equation of the state of the con- 
fined fluid. Classical density functional theory [22] 
can also be applied to study the slit pore geometry, 
with very good agreement with experiments and 
simulations. Though the theoretical work devel- 
oped in these pages is not a competitor of density 
functional theory, it has the advantages of having 
analytical expressions, and the possibility of eas- 
ily introducing two immiscible components: for in- 
stance one or two layer lubricants wetting the walls 
and a gas or a liquid filling the rest of layers forming 
the capillary volume. 
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Appendix A 

Expressions of quantities used in Eqs. 5-8: 
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